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Abstract Recent experimental work has shown that eukaryotic cells can swim in a fluid as well as crawl 
on a substrate. We investigate the swimming behavior of Dictyostelium discoideum amoebae who swim by 
initiating traveling protrusions at the front that propagate rearward. In our model we prescribe the velocity 
at the surface of the swimming cell, and use techniques of complex analysis to develop 2D models that 
enable us to study the fluid-cell interaction. Shapes that approximate the protrusions used by Dictyostelium 
discoideum can be generated via the Schwarz-Christoffel transformation, and the boundary-value problem 
that results for swimmers in the Stokes flow regime is then reduced to an integral equation on the boundary 
of the unit disk. We analyze the swimming characteristics of several varieties of swimming Dictyostelium 
discoideum amoebae, and discuss how the slenderness of the cell body and the shapes of the protrusion 
effect the swimming of these cells. The results may provide guidance in designing low Reynolds number 
swimming models. 

Keywords Low Reynolds number swimming • self-propulsion • amoeboid swimmimg • metastasis • robotic 
swimmers 
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1 Introduction 


Cell locomotion is essential throughout the development and adult forms of uni- and multi-cellular organ¬ 
isms. It is beneficial in various types of taxis, in morphogenetic movements during development, and in the 
immune response and wound healing in adults, but it plays a deleterious role in cancer metastasis. When 
movement is in response to extracellular signals, it involves the detection and transduction of those signals, 
which can be biochemical, mechanical or of other types, the integration of the signals into an intracellular 
signal, and the spatio-temporal control of the intracellular biochemical and mechanical responses that lead 
to force generation, morphological changes and directed movement ( Sheetz et al.\\999 \. Many single-celled 
organisms use flagella or cilia to swim, and many mathematical models of swimming in such organisms have 
been developed ( |Lau g a & Powers| | 2009[ ). The movements of eukaryotic cells that lack such structures fall 
into two broad categories: mesenchymal and amoeboid ( Biname et al.|201^ . The former can be character¬ 
ized as ‘crawling’ in fibroblasts or ‘gliding’ in keratocytes, and involves the extension of either pseudopodia 
and/or lamellipodia driven by actin polymerization at the leading edge. This mode dominates in cells such 
as fibroblasts when moving on a 2D substrate. On flat surfaces the predominant protrusions are lamellipodia 
and these processes suffice, but in the extracellular matrix (ECM) the protrusions and cell body are more 
rounded, and the cells may also secrete matrix-degrading proteases (MMPs) and ‘tunnel’ their way through 
the ECM ( Martins & Kolega|2006 Mantzaris et q/.|2004[ ). 

In the amoeboid mode cells are more rounded, and move through the ECM by avoiding obstacles when 
possible, rather than removing them. Thus they rely less on attachment to the ECM and degradation of it, 
and instead exploit variations in the ECM to move through it by shape changes. In this mode force trans¬ 
mission to the extracellular matrix (ECM) depends on shape changes driven by localized remodeling of 
the cytoskeleton (CSK) and myosin contraction ( |Insall & Machesky|2009| l. Cells such as leukocytes, which 
normally use the mesenchymal mode in the ECM, can migrate in vivo in the absence of integrins, using a 
’flowing and squeezing’ mechanism ( Lammermann et al.|2008| l. The human parasite Entamoeba histolytica 
uses an extreme form of amoeboid movement called blebbing, in which the membrane detaches from the 
CSK (Maugis et al. 2010| l, whereas zebraflsh primordial germ cells and Dictyostelium discoideum (Dd) cells 
move using a combination of blebbing and protrusion by spatio-temporal control of the membrane attach¬ 
ment to the CSK (Blaser et al.||2006 |Diz-Munoz et al.||2010 Yoshida & Soldati]|2006[ Zatulovskiy et al. 


20141. Cells using this mode can move up to forty times faster than those using strong adhesion (Renkawitz 


& Sixt|2010|l. Recent experiments have shown that numerous cell types display enormous plasticity in lo¬ 


comotion in that they sense the mechanical properties of their environment and adjust the balance between 
the modes by altering the balance between parallel signal transduction pathways ([Renkawitz et a/.||2009 


[Renkawitz & Sixt]|2010| l. Thus crawling and swimming are the extremes on a continuum of locomotion 
strategies, but many cells sense their environment and use the most efficient strategy in a given context. 

Heretofore mathematical modeling has focused primarily on either the mesenchymal mode, in which 
cells crawl via attachments to a substrate (Danuser et a/.|2013]l, or on microorganisms that swim using flag 


ella or cilia ( Suarez & Pacey|2006[ Berg & Anderson|1973 Lowe et a/. [1987 Gibbons|1981 Sleigh et al. 
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1988 Ishimoto & Gaffne^|2014|l. Here we analyze swimming of larger cells, motivated by recent exper¬ 


iments which show that both neutrophils and Dd can swim - in the strict sense of propelling themselves 
through a fluid using only fluid-cell interactions - in response to chemotactic gradients ( |Barry & Bretscher] 
2010[ Bae & Bodenschatz|2010 [Van Haastert|20ir| l. A basic question that arises is what pattern of shape 


changes are effective in propelling a cell. It has been reported ( Salbreux et a/.|2007) l that freely suspended 
Swiss-3T3 fibroblast cells without cell-substrate adhesion can exhibit oscillatory shape dynamics, and some 
of them may also exhibit periodic bleb dynamics correlating with the oscillations, where blebs are hemi¬ 
spherical membrane blisters induced by cortical contraction ( |Fackler & Grosse|2008| |Paluch et a/.||2005] l. 
The combination of both dynamics may result in a random oscillation of the cell in space. In comparison 
with such random motion, Dictyostelium and neutrophils can utilize the amoeboid swimming mode in which 
the cell body is elongated and small protrusions that provide the momentum transfer needed for motion are 
propagated from front to rear ( Barry & Bretscher|2010[ Bae & Bodenschatz|2010[ Van Haastert|201 1) . 

In the following section we formulate the basic problem of swimming by shape deformations in 2D at 
low Reynolds number, which is the relevant regime for single cell movement. In 2D one can introduce a 
stream function, which leads to a biharmonic equation, and the general solution of the Stokes problem is 
expressed in terms of two analytic functions - the Goursat functions - that are determined by the motion 


of the boundary of the swimmer (Bouffanais et al. 2013 Chambrion & Munnier 2011 Gherman et al. 


20001. This in turn leads to an integral equation for one of these functions, and the second function can 


then be expressed in terms of the first. In Section]^ we study the motion of Dd in 2D, and approximate the 
shape changes using polygonal approximations. Using the Schwarz-Christoffel transformation, we reduce 
the problem to the solution of a linear system of equations for basis functions on the boundary of the unit 
disk. We show that realistic propagating shapes can produce propulsion at speeds in the range observed 
experimentally using realistic choices of the parameters. 

Since cell movement, whether on a solid substrate, in a fluid, or in a complex medium such as the ECM, 
involves the interplay between biochemical and mechanical processes, dissecting the roles of each is a first 
step toward an integrated description of movement. A major objective of our work is to answer a question 
posed by experimentalists, which is ‘How does deformation of the cell body translate into locomotion?’ 
( Renkawitz & Sixt|2010| l. A longer-range goal is to produce a unified description for swimming that inte¬ 
grates signaling and mechanics in viscous and viscoelastic environments similar to the extracellular tissue 
environment. As one experimentalist stated ’the complexity of cell motility and its regulation, combined 
with our increasing molecular insight into mechanisms, cries out for a more inclusive and holistic approach. 


using systems biology or computational modeling, to connect the pathways to overall cell behavior’ (Insall 
& Machesky|2009 1. 


2 Low Reynolds Number Swimming Problems 

2.1 General description of swimming mechanics at low Reynolds numbers 


Recent interest in the motion of biological organisms in a viscous fluid was re-kindled by Purcell’s 1977 
description of life at low Reynolds number (LRN) ( Purcell|1977 1, and a wide variety of applications have 
been analyzed since then. A review of some of these is given elsewhere ( |Lauga & Powers||2009| l, and we 
only describe the relevant background. We consider motion in an incompressible Newtonian fluid of density 
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p, viscosity p, and velocity m, for which the governing equations are 

du 

+ p(«- V)m = V-T + /= -Vp + ^Au + f, 
at 

Vm = 0 


( 2 . 1 ) 

( 2 . 2 ) 


where T = —p5 + p(Vm + (Vm)^) is the viscous stress tensor, f is the external force field, which we assume 
to be zero or include it in the pressure hereafter. This is appropriate for the experimental configuration 
described in ( |Barry & Bretscher|20101 l, since the fluid had a variable density that allowed cells to And a point 
of neutral buoyancy. The Reynolds number based on a characteristic length scale L and speed scale U is Re 
= pLU /p, and when converted to dimensionless form and the symbols re-defined, the equations read 


du 

ReSl^ — \-Re(u-'V)u = —Vp + Au 
at 

V.M = 0. 


(2.3) 


Here SI = (OL/U is the Strouhal number and O) is a characteristic frequency of the shape changes. When 


Re ^ 1 the convective momentum term in (2.3 i can be neglected, but the time variation requires that Re SI = 


2013 


Wang & Ardekani 20121, but 


(OL? jv 1. This is not always true, even for small swimmers (Ishimoto 
as we show next, it applies here. 

The small size and slow speed of cells considered here leads to LRN flows, and in this regime cells 
move by exploiting the viscous resistance of the fluid. For example, Dd amoebae have a typical length 
L ~ 25pm and can swim at U 
lO^kg m“ 


3pm/min (Van Haastert 20111. Assuming the medium is water (p 
“^Pa-s), and the deformation frequency o ~ 1/s, Re ~ €?(10^®) and SI ~ 


p ~ 10^^Pa-s), and the deformation frequency o ~ 1/s, Re ~ and SI ~ In 

fact the experiments are done in oil that is significantly more viscous ( Barry & Bretscher|2010| l. In any case, 
when both inertial terms on the left hand side of ( |2.3| l are neglected, which we assume hereafter, the flow is 
governed by the Stokes equations 


pAM-Vp = 0, 


V.M = 0. 


(2.4) 


Throughout we consider the propulsion problem in an infinite domain and assume that the fluid is at rest at 
infinity. 

Since time does not appear in the equations, a time-reversible stroke produces no net motion, which is 
the content of the famous ‘scallop theorem’ ( Purcell|1977 1. Furthermore, since we assume that there is no 
inertia in the fluid (momentum transfer is assumed to be instantaneous), in the Stokes regime there is no net 
force or torque on a self-propelled swimmer moving in an infinite fluid that is at rest at infinity. One can see 
this by integrating the stress over the boundary of the swimmer and then equating this via the momentum 
equation to the stress on a circle at infinity ( Shapere & Wilczek||1989b l. Therefore movement is a purely 
geometric process: the net displacement of a swimmer during a stroke is independent of the rate at which 
the stroke is executed when Re = 0, and approximately so as long as Re SI remains small enough. 

Let i2{t) G R" be the swimmer (a compact set with a sufficiently smooth boundary). Throughout we 
consider a fixed global reference frame (x,f) and a body frame {X,t) attached to the swimmer, where time 
is scaled the same in both frames. We use the notation dQx{t) and dQx{t) to denote the boundary of £2 
in the fixed and body frames, resp., and when observed from the body frame, the shape deformations of 
the swimmer are specified as Us{X,t) for X G dQx{t). Denote the uniform (rigid-body) translational and 
rotational velocities of the swimmer when observed from the fixed frame as U{t) and respectively. 
Then the instantaneous velocity on the swimmer’s surface dQx{t) is 

u{xp) = U {t) + 0}{t) X X + Us{X ,t) 


( 2 . 5 ) 
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where u satisfies ( |2.4| l. 

A swimming stroke is specified by a time-dependent sequence of shapes, and it is cyclic if the initial 
and final shapes are identical ( |Shapere & Wilczek|1989b] l. The canonical LRN self-propulsion problem is 
- given a cyclic shape deformation specified by Us, solve the Stokes equations subject to the zero force and 
torque conditions 

F{t) = f T-n=0, T{t) = f jeA(T-n)=0 (2.6) 

JdQx{t) 

and the boundary conditions 

u\dQ^Us + U -ymxx, u\x^c<,^0 (2.7) 

where n is the exterior normal. 

How does the swimmer move if force- and torque-free? It would appear that at best it can only translate 
or rotate at a constant rate. However, an arbitrary deformation of the surface will not satisfy p.6| l in general, 
and in particular, will generate a flow at infinity. However, since we assume that the fluid is at rest there this 
flow must be counteracted by an imposed flow, which then defines the rigid linear and angular velocities 
of the swimmer, and leads to satisfaction of the force-free and torque-free conditions (| Shapere & Wilczek] 
1989b] l. Since the shape changes are time-dependent the counterflow is also, as are the rigid translations and 
rotations. 

Thus the canonical LRN swimming problem is: given a cyclic sequence of shape deformations by speci¬ 
fying Us on the boundary d£lx of the swimmer, solve the Stokes equations \2.A\ for a trial velocity field u and 
then use the force-free and torque-free conditions to determine U and CO so as to satisfy the boundary con¬ 
ditions at infinity. In what follows we restrict attention to swimming in two space dimensions, and describe 
the problem as the 2D LRN swimming (2DLRNS) problem. 


2.2 The 2D swimmer 


In two space dimensions techniques from complex analysis can be used to significantly simplify the problem 
of computing solutions to the Stokes equations for LRN swimming problems. Two main methods have 
been developed that can be applied, one that was first developed by Muskhelishvili to solve problems in 
elasticity (|Muskhelishvili| 19'771, and another that is essentially a boundary integral method! Pozrikidis|1992[) 


for 2D problems ( Greengard et a/.|1996[ Kropinski|1999[ 200 1[ 2002| Kropinski & Lushi|201 1 1. Significant 
analytical insights can be gained using Muskhelishvili’s method, including the application of control theory 
to LRN swimming (Shapere & Wilczek 1989b|a| Kelly 1998[ Kelly & Murray|2000) l, and we use this method 
in this paper. 

In 2D the incompressibility condition V •« = 0 in (2.4 1 can be satisfied by introducing a stream function 
A (z,z;f), which is a real-valued scalar potential such that 

dA ,dA 
dy dx 


Here and hereafter we use m S C to denote the velocity field (denoted by u in Section 2.1 1 in the complex 
z-plane. Then the Stokes equations \2.A\ imply that A satisfies the biharmonic equation 


= 0 . 


( 2 . 8 ) 


The general solution of ( |2.8| l can be expressed by Goursat’s formula ( Muskhelishvili|1977 1 

A{z,z;t) =3i[z(j)iz;t) -\-xiz;t)] 


(2.9) 
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where for any t, (l>(z;t) and x(z;f) are analytic functions on the fluid domain C\i2(t) and continuous on 
C\intf2(f), where intf2 denotes the interior of and x(z;f) are the Goursatfunctions. To simplify 

the expression of physical quantities and the following discussion, we hereafter impose the substitution 
(j) —> —icj), X suggested by Shapere and Wilczek ( Shapere & Wilczek|1989b| l, and Table[^gives 

the expressions of several physical quantities in terms of these functions. In the table and hereafter we denote 
by ' for simplicity, n = —idz/ds gives the exterior normal to dQ {i.e., directed into the fluid domain), and 
s denotes the arc length, traversed counterclockwise. 


Table 1 Representation of various physical quantities by the Goursat functions. 


Velocity 

u = 0(z)-z0'(z)-l'(z) 

Pressure 


Vorticity 

^ = -43{0'(z)} 

Stress force 

/ = 4fi3i{(l)')n - 2ii{z^" + X")n 

Stress force 
(differential form) 

fds = -2iiid (0 + Z0' + X') 


For swimming problems in a 2D Stokes flow in the unbounded domain C/Q{t), we require that the 


stress vanish at infinity, and as a result the Goursat functions must take the general form (Muskhelishvili 
1977[ Greengard et q7|199^ : 


0(z,f) = - 


X{t) + iY{t) 


27r(l + K-(f)) 


logz + 0(z,t) 


¥iz,t) = x'(z,t) = 


27r(l + K-(f)) 


logz + v/(z;f) 


( 2 . 10 ) 


( 2 . 11 ) 


where ((>(z,t) \(f(z, t) are single-valued and analytic on C /Q (where C = C U1°°}). For the self-propulsion 

problem we must also require that X = Y = 0 to ensure a bounded velocity at infinity. We then compute the 
translational and rotational velocity for a trial pair (0,x), and when these are subtracted from the flow the 
motion is force-free and torque-free ( |Muskhelishvili|1977[ Greengard et a7|1996| l. Thus for 2D swimming 
problems, the Goursat functions 0(z) and v/(z) should be single-valued and analytic on C/f2, and therefore 
they have Laurent expansions in C /Q of the following form. 


0(z,f) = ao{t)- 

Wizd) = boit)- 


a-2{t) 

z zf 

b-\(t) b_2{t) 


( 2 . 12 ) 


z 


(2.13) 
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As a result, the 2DLRNS problem with specified shape changes as described above can be rephrased as 
follows - find functions <p(zf) and that are analytic on C/Q{t) and continuous on C/int Q{t) such 

that 


<j){z,t)-z<j)'iz;t)-\if{z,t) = V{z,z,t) {zGdQ{t)). (2.14) 

Here V{z,z,t) for z G dQ(t) is the velocity boundary condition determined by the shape changes, namely. 


the complex form of Ms as introduced in Section 2.1 Equation (2.141 will be referred to as the boundary 


condition constraint on (j) and y/ in the study of 2D LRN swimming. 


2.3 Pull-back of the problem to the disk and derivation of the Fredholm integral equation 


As we mentioned earlier, there are two main methods for solving the 2D LRN problem - which is to say to 
solve ( 2.14| i - Muskhelishvili’s method and the integral equation method. In the integral equation method 
the velocity boundary condition leads to an integral operator defined on the swimmer’s boundary, while in 
Muskhelishvili’s method the approach is to first map the z-plane to a fixed complex computational i^-plane, 
on which the integral operator is applied. In the first case one deals with a moving boundary whose stress 
field is specified, while in Muskhelishvili’s method we consider velocity boundary condition and the problem 
can be pulled back into a fixed boundary problem. Generally speaking, the integral equation method is useful 
when the stress field along the boundary is known or multiple bodies are involved. On the other hand, for 
a single deformable swimmer it is easier to treat a sequence of shape deformations using Muskhelishvili’s 
method, which also facilitates the use of control theory to 2D LRN swimming systems and simplifies the 
design and study of micro aquatic robots. Hereafter we focus on the use of Muskhelishvili’s method. 

Suppose that the cell occupies an open, simply-connected, bounded region i2(f) in the complex z-plane 
at time t. Let D = {i^GC:|i^|<l}be the unit disk in the computational i^-plane. The regions exterior 
to D and to Q{t) in the extended complex planes, i.e., C/D and C/ Q(t), are both simply-connected as 
the infinity point o° is included in C. The Riemann mapping theorem ( Ahlfors|l978 1 ensures the existence 
of a single-valued analytic conformal mapping z = w(^;t) which maps C/D one-to-one and onto C/Q{t), 
and preserves the cor respondence o f infinity, i.e., w{o°-,t) = Moreover, this mapping can be extended 
continuously to C/D (Younes 20101, and w maps d{C/D) = 5' to dQ{t). 

We assume that the mapping preserves the point at infinity, and then have the following result regarding 
the form of its Laurent expansion ( Ahlfors|1978| l. 

Lemma 1 Suppose that Q{t) is a non-empty open bounded simply-connected domain in C, and that z = 
w(^;f) is a conformal mapping from the exterior of the unit disk D to that preserves the point at infinity. 
Then w{^',t) has a Laurent expansion of the form 


X XX «-l(0 CC-lit) a-nit) 


(2.15) 


C ' C" 

where ai{t) ^ 0 and |^| > 1. 

When we consider in the body frame, Q (f) gives the configuration of the swimmer, in which case we require 
Oo = 0 and a^i G . 


in (2.15 1 . 


Hereafter is the conformal map which maps C/D onto C/Q{t) such that w(oo;f) = extended 

to the boundary of the swimmer, which is given by dQ (f) = {z(f) = w((7;f); O’ G 5*}. We impose a no-slip 
condition on the boundary, and therefore 




















Qixuan Wang, Hans G. Othmer 



Fig. 2.1 The conformal mapping from C/D to C/f2. 


Suppose that 0(z,t) and are the solutions to the boundary condition constraint given by (2.14i. 


Let 


0(C;f) = 0WC;O;O, TO = r(w(C;0;0 (ICI>i) 


thus and are functions on -plane which are analytic on C/D and continuous on C/D at any 

time t. Hence on \ Q ^ 1, we have Laurent expansion for <l> and ^ given by 


0(C;O + + + 

(2.17) 

TO0=fio(0 + ^ + ^ + --- 

(2.18) 

where (A: = 0, — 1,—2, •••) are continuous functions from [0,°°] to C. 

If we substitute these expressions into (2.14 1 and omit the time index f, then the velocity field can be 
expressed as a function on the i^-plane as 

«(C, C) := «(w(C),TO) = ^(C) - - TO 

w'(C) 

(2.19) 

where ' represents d;^. If we let y(c7;t) denote the boundary condition of the velocity in 
boundary condition constraint (|2.14|) can be generalized to the condition 

general, then the 

w'(cr) 

(2.20) 


In the 2DLRNS problem the boundary condition is given as in ( |2.16| l, i.e., V{o\t) = (9,w((7;f), but in the 
discussion below we will adopt the more general form V (<7;f) = Ln 
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Equation p.20| l involves both <l> and but this can be reduced to an integral equation in (p alone 
using a technique developed in Muskhelishvili (19771, where it was originally developed for the solution 
of problems in elasticity rather than the Stokes flows, and the resulting integral equation is the same. We 
present the precise statement here for completeness. The reduction relies on a decomposition of the boundary 


condition y(<7) derived from the following Plemelj formula (Sokhotskii 1873 Cima & Ross 2006 England 

[ 2012 ] ). 


Theorem 1 (Plemelj Formula) Suppose V (cr) is continuous on = {cr € C; |(7| = 1} and, for a particular 
Ob S S\ satisfies the Holder condition 


|y(o)-y(oo)| <c|o-oo|“, oe5^ 


( 2 . 21 ) 


for some positive constants C and <X. Let 


v{Q 


f v{<y) 
2ni Js' o - C 


da 


( 2 . 22 ) 


for ^ G C/5'. Then the limits 

y^(oo):= lim y(roo) and y^(oo) := lim y(oo/r) (2.23) 

r -!-1 “ r -!-1 “ 

exist and moreover, y ^ (oq) — y^ (Oo) = V (Oo). Furthermore, the Cauchy principal-value integral 

RV. I := lim f 

Js^ (y — CTq (7Q|>e (7 — do 


exists and 


y+(oo)+y (oo) 


1/..V. f 

Tti Js' o — Oo 


We as sume that the boundary condition y(o) in the swimmin g pro blems satisfies the Holder condition 
as in ( 2.21 1 for any Oo G 5*. Eor the function V{Q defined in (2.22i, let = V{Q for |i^| > 1, and 

V^{Q = y(C) for ICI < 1’ then V^{Q is analytic for |i^| < 1, while is analytic for \Q > 1 with 

y+(Q —0 as —» °o. M oreover, according to the Plemelj formula, both y+ and can be continuously 
extended to 5' by (2.23 i, and we have the decomposition V = on 5'. 


T o continue, we hrst intr oduce some notation. Eor an analytic function f{Q, we dehne /(C) as /(C) = 

/(C) (Muskhelishvili 1977 i. With this notation, suppose that /(C) is analytic on |C| ^ R for some R > 0, 
then /(I / C) is analytic on |C| The proof of this is straightforward. Without loss of generality we assume 
that /(C) is analytic on |C| > R, and then on |C| > R we have the following Laurent expansion for /(C) 

/(C)=/o + ^ + ^ + --- + ^ + --- 


/(^) —/o+/lC+/2C^ H-f/nC"H- 


Hence 
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which is easily seen to be analytic on 1| < 7?. By definition of the function /, we have 


/(^)—/(=)—/o+/iC+/2C H-l■/«C H-—/o+/iC+/2C^H-+ ■ 

which is easily seen to be analytic on 1 1 ^ | < R as well. This establishes the assertion. 

To derive an integral equation for ^>, let | ^ | > 1 and apply the functional operator 


— / — 

Ini Js' <7- C 


da 


to both sides of (2.20i. is analytic on |i^| > 1, and continuous on |i^| > 1, so 


<P{a] 


da = -^(C) 


2ni iji a - C 

for Id > 1. Moreover, 'f^( 1 / d is analytic on | d < 1 and continuous on | d < 1, and thus 


^{a) 


da 


-da — 0 


2ni Js' cj - C 2ni (y-C 
On the other hand, for the right-hand side of ( |2.20| l by Plemelj formula we have 

2niJsia-i; 

Thus the final result is Fredholm integral equation 


^(0 + ^- 


w{a) <P'{a) 


2niJs^ w’(a) c^-C 


da = -V+{Q (ICI>1) 


(2.24) 


where —y+(d is the analytic part of V(d when |d > 1, namely, —y+((j;f) = Y.n<o^n{t)a'‘. This equation 
has only one unknown function <l> and once it is known "F can be obtained from ( 2.20|l. Th e joint solution is 
unique in the sense that the constant terms Aq and Bq of <P and "F in equations ( 2.17|2.18 1 may vary, but the 
difference Aq — Bq is uniquely determined. 


2.4 Expressions of physical quantities by the pull-back of Goursat functions 


Given the pull-back of the Goursat functions determined by (2.24), whose Laurent expansions are in the 


forms given in equations (|2.17|2.18]l, we can obtain the expressions of several physical quantities of interest 


in a typical swimming problem (Shapere & Wilczek 1989b Gherman et al. 2000[ Shapere & Wilczek|1989: 


Avron et q/.|2004| i. In the following discussion we scale the length by R and the time by T, where R usually 
corresponds to the radius of the cell when it is in the shape of a disk, /.e., nR?' = Area of the cell. 


1. Rigid motions (I): translation. Following the approach ussed in ( Shapere & Wilczek|1989b| l, let Uoc and 
Woo be the translational and rotational components of the far field behavior when prescribing boundary 
condition ( 2.16| l, respectively; then the actual translational velocity U and the rotational velocity (O of 
the swimmer are given by 


U = -t/ 00 , CO = -Woo 


(2.25) 
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as the fluid is static at far field. As shown in appendix^ t/«, can be computed from the relation 

Uoo=ao — bo=AQ — Bo (2.26) 

where the ag, bo are the leading order terms of (j) and xf/ as given in equations ( 2.12|2.13 1 , and Aq, Bq 
are the leading order terms of <l> and f' as given in equations ( 2.17|2.18 i. Then the net translation of the 
swimmer at any instant t is given by 


Tr{t)= f -U^{t)dt= f \-Ao{t)+Bo{t) 
Jo Jo L 


dt 


and therefore the average velocity of the swimmer within a period isU = Tr{T}/T. 

2. Rigid motions (II): rotation. The rotational velocity co can be obtained from the following relations. 


CO = 


T{V;w) _ -47r/r3/7_i _ -47r/r3(B-iai) 


7’(y™‘;w) T{V’^°\w) 




(2.27) 


Here b-\ and B \ are the coefficients of and resp., in the Laurent expansions of i/(z) and 
in (2.13 I and (2.181, resp., and is the coefficient of the leading order ^ term in the conformal 
mapping z = vr(Q in (2.15 i. Further, T{V',w) is the torque resulting from the current shape z = vr(Q 
and deformation y((T) of the swimmer, while ^(y™*; w) is the torque resulted from a rigid rotation of 
a swimmer also with shape z = vv(i^) but a uniform rotational velocity field y™'^((7;f) = /w((7;f). More 
detailed analyses of U and co are given appendix^and appendix]^ 

3. Force distribution. From Table[^ we see that fds = —2ijid{2^ — y) on dQ. The pull-back of the force 
distribution is: 




(2.28) 


4. Power expenditure. The power expenditure is calculated by integrating the stress times the velocity on 
the surface of the swimmer: 


^ = ufds = -2ii^ [ V(a)(2<P'(a)-V'(a)) 

Jdn Js' 


da 


With boundary condition generally given as V (a) = together with the expansion of ^{a) as 

given in p.l7|i, we have 


' = Anp L « 

n>\ 




We define the average power expenditure within a period to be 

Ip [ ^it)dt 

I Jo 

5. Performance. We define the performance of the swimmer as 


Tr{T) 

^ ll^iiidt 


E=L = 


(2.29) 


(2.30) 


(2.31) 


E measures the distance traveled in one period divided by power expended in a period. 
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6. Area of the swimmer. We usually require that the total mass of the swimmer be constant for cells swim¬ 
ming through a fluid by shape changes. For a swimmer of constant density in 2D this becomes an area 
conservation constraint. Suppose that the shape changes of the swimmer are given by ( |2.15[ ), then the 
swimmer’s area is 

Area(f) = ^3 ^ wdw. 

By a direct calculation we obtain 

Area(f) = 7r(|ai|^- |a_i|^-2|a_2|^-- ) (2.32) 

For incompressible swimmers, we require that Area(f) = Constant. 

A list of scales and units of these physical quantities is given in Table]^ 


Table 2 A list of physical quantities 



Notation 

Scale 

Unit 

Length 


R 

pLm 

Area 


Rf 

llm^ 

Time 

t 

T 

s 

Mean velocity 

U 

R/T 

IXm/s 

Force density 

f 

M/r 

pN / pLm 

Mean power 


PlR^/T^ 

pN ■ jXm/s 

Performance 

E 

T/W 

l/pN 


3 Shapes for conformal mappings with finitely many terms 


In general the Fredholm integral ( 2.24| i cannot be solved analytically for an arbitrary conformal mapping 
w{Q. However when w{Q has only hnitely many terms, we hnd that (2.24i can be analytically simplihed to 
a linear relation between the coefficients of 0 and F+. In view of the fact that any conformal mapping can 
be approximated by truncating its Laurent expansion, this approach is sufficient in the study of 2D Stokes 
flow swimming problems. 

We consider a sequence of shape changes whose corresponding conformal mappings always have Lau¬ 
rent expansions up to (—A)-th order. Then 


^ /N OC-i(t) a-2{t) OC-N{t) 

w(C;f) — ai(f)C-f ao(f)-l- - - \ - ^ -1-1-^— (3.1) 

In general, the term with n > 0 gives n + l angles along the periphery of the swimmer. Figure [ lT] gives 
an illustration of the relation between the conformal mapping and the shape of the swimmer. 

We assume that the boundary condition is in the general form of V((7;f) = LA„(f)(7", and we hrst prove 
the following result in some particular cases. 
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Fig. 3.1 The relationship between the conformal mapping and the shape of the swimmer: a): Shapes determined by w{Q = 4iJ + 
where N = \ (solid line, ellipse), N = 2 (dashed line, triangle), or A = 3 (dotted line, quadrilateral), b): A sequence of shape changes 
defined by = 3^ + cos(2m)if — sin(2m)^“^. 


Lemma 2 Suppose that the conformal mapping w{l^) is given by ( |3.1| ) and the boundary condition is given 
as V (cj) = Let —y+(^) be the analytic part ofV (^) when |^| > 1, which has the form —y+((J;f) = 

£„<q A„(f)(7". If there exists a function /iv(C) analytic on |^| < 1, continuous o« |^| < 1, and 


fN{<y) 


CTGS* w'((7) 


for some N € N >2, then for any —y+ (Q with Laurent series such that A_i — A_2 — • • • — — 0 

for N > 2 (ifN = 2 there are no restrictions on the coefficients), the solution to the boundary value problem 
(|2.20|) is given by 


0{Q = -y+(C) 

W (C, j 


for\C\ > 1- 

Lemmaj^can be proved by direct calculation as follows. 
Proof The function —y+ has the representation 


-^+(0 


■^-(A-1) 

I^A-l 




+ ••• 


Suppose that ^’(C) = Then on the unit circle 5', we have 

^>'(a) = -a^[{N- 1)A_(2^,_i) +A^A_At(7H-l-nA_„c7"^^+' H—] 

and the integral term in ( |2.24| l becomes 

J_ r 

27tiJs^ w'{a) (T-C 

^ ^)'^-(A-i)-l-nA_„c7"^^+* H— ]da. 

Since the function 


— fN{^)[{N — l)A_(^_i) +A^A_jv(JH-|-nA_„(j" H-] 
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is analytic on |ff| < 1 and continuous on |(j| < 1, and since \ Q > 1, the integral 


^ f g(g) 

27Zi Js' cr - C 


da = 0 


by the Cauchy Integral Theorem. Hence the left hand side of (2.24 1 reduces to ^’(C) = substi¬ 

tuting this result into (2.20i, we obtain the expression of 'F as in the lemma. Q.E.D. 


With w given by (3.1 1 we have 


w{a) _ aia + oco + a-ia ' + a-2(y - ha-^a 

w'{a) ~c^-al\a^-2a~2(y^ - 




hence we may take 


/tv(C) = c 


^aiC + cto + ct-iC ' + 0!_2C -^ 


ai -a_iC^-2a_2C^--^o:_7vC^^' 

a-N + -f • • • + * -f 0!oC^ + Ctl 


«! - -2a_2C^- Nu-nC 


N+l 


(3.2) 


If no singularity of fy lies inside the unit disk |^| < 1, then is analytic on |i^| < 1 and we may apply 
Lemma|^to obtain the solution of the boundary condition constraint ( |2.14| l. 

The Laurent expansion of any function ‘P{C) that is analytic on |i^| > 1, continuous on |i^| = 1, has the 
form of ( |2.17[ ). Let be the set 




{^(C) 


C C' 


1-3 


^3 


+ ■■■ ;l^^|A-i| < +°°| 

k>i ^ 


of all functions ‘PiQ whose coefficients of the Laurent expansion satisfy Y,k>i ^|A-/i| < +°° and vanish at 
inhnity. is a Banach space endowed with the norm || <l> \\ = Lk>i k|A_i| (IChambrion & Munnier|2011|. 


Now with w given by (3.1 1 , we dehne a conjugate linear operator K{w) on ^ as 

w(( 7 ) <P'{a) 


[Kiw)o<p]iC) = — 


27 ziJsi vv'(a) 


da 


(3.3) 


Let = {Aq G C}, then ^r(w) can be extended to such \h&tK{w)\y^ = 0 in accordance with (3.3 1. 

Now for any e ^ 0 (2.24 i can be written as 

[I + K{w)\^ = -V+ 

where I is the identity operator. Lor n e Z+, let 




and 


It is easily seen that ^„^’s and ^„^’s are real 1-dimensional linear subspaces of 5^. With /^v given by (3.2i, 
by Lemma]^ ^r(w) vanishes on each and for n =N— 1,N,N+ 1, • • •, and thus each 5^^ or is 
invariant under the operator I + .K’(w), i.e., 


-K{w)] 


-^r(w)] 1^3 — I yS 
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Define 


-(N-l) 


Span^ji^ 


n=— 1 


which is a real 2{N — 2)-dimensional linear subspace of and 

consists a basis of If /n is analytic on |i^| < 1, then the action of K{w) on this basis can be expressed as 

fN{<y) 1 _ n 




's' (y-C 

fN{(y) 1 


da = 


n 


_ ^ Jn I'k 


da = —■ 


k\ 


( 0 ) 


k=0 

in ^ 


'N-n- 


I L 


k\ 


l-k 


2nJs' a-C 

for n = 1,2, • • • ,— 2. Notice that the sum in the above equations is simply the hrst N — n — I terms of 
the Laurent expansion of /a'(C) at = 0. It is easily seen that on =5^, the operator K{w) has a matrix 
representation as 

[^:(w)] o (C-', C-^ • • •, = K{w) ■ 

where K(w) is a {N — 2) x {N — 2) matrix with the form 


( fr^\o) 

{N-3)\ 


K{w) = 


(N-4)! 


(N-4)! (N-5)! 






(N-5)! (N-6)! 


fm MY mA 

2 ! 1 ! 0 ! 

2 /;( 0 ) 2fy{Q) 


V 


(jV-3)/;(0) (jV-3)/^(0) 
(iV-2)A(0) 


0 ! 


0 ! 

0 


1 ! 

3/v(0) 

0 ! 

0 


0 

0 


0 ! 

0 


0 

0 


0 

0 


(3.4) 


This can be summarized as follows. 


Lemma 3 Suppose that the conformal mapping w(Q is given by (3.1 i and the boundary condition is given 
as V(<7) = If there exists a function /iv(C) analytic on |^| < 1, continuous on\l^\< 1, and 


fN{<y) 


aeS' w'{a) 


for some N G N >2, then the integral ( 2.24fl reduces to linear relations between coefficients of and 
V: 
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1. For n = 0 or n>N — = A_„; 

2. For 1 < n < N — 2: 


{In-2+K{w)) • • ,5RA_(^_2))^ = 

(4^-2 —(3A_1, • • • ,3A_(^_2)) = !3AAr_2) 


where In- 2 is the [N — 2)-identity matrix, and K{w) is the (N 


2) X (N — 2) matrix given by {3.4 


Lemma provides an algorithm for solving the Stokes equation of an infinite 2D Stokes flow when 
the shape deformations have finitely many terms in the conformal mapping w{Q. Once we solve for <l> 
according to Lemma|^ ^ is given by the following expression. 


«! 

I 

«r 


= -^ + (A_i-a_i)C + (A_2-o:-2)C^H- h (a_a, - a_A-)C^ 


^ + a_iC + a-2C" + --- + «-ivC'^ . 2A_2 A^A_^ 


«! - 


a_i 2a_2 


^2 


^3 


Na-N 
■ ~^N+r 


^3 


|XA'+1 


(3.5) 


4 Modeling of swimming Dictyostelium amoebae 


As mentioned in the Introduction, Dictyostelium amoebae can move in a fluid environment by a combina¬ 
tion of blebbing and protrusions, either of which involve rapid shape changes and neither of which require 
attachment to a substrate. Experimental observations on the movement of Dd cells reported in [Barry & | 
Bretscher ( 2010| l, Bae & Bodenschatz (20101 and Van Haastert (20111 have recorded cell shape changes, 
speeds, and periods of the cyclic motion, and we use their data here to compare with theoretical predic¬ 
tions. The movement usually involves protrusions that are initiated at the leading edge of the cell and which 
propagate toward the rear. Typically the cell body is elongated, and multiple protrusions propagate along the 
cell, van Haastert ( Van Haastert|201 Ij l reported an average of three protrusions, as illustrated by the cartoon 
model in Figure 4.1 (a), while from the experimental images (Figure |4~T| (b)) of a swimming Dictyostelium 
reported in Barry et.al. ( [Barry & Bretscher|2010| l, we see that one protrusion travels along one side of the 
cell and disappears at the rear of the cell, then another protrusion appears on the other side and repeats the 
process. 

Van Haastert ( [Van Haastert|201 1) observed that the protrusions travel directly down the cell body and not 
in a helical fashion. Thus there is no clear evidence that the cell is rotating around its symmetry axis, and as 
a result we consider the 2D model developed in the preceding sections to be a reasonable simplification of a 
3D swimming cell. In Section |4T[ the numerical scheme used to constmct the shape of the cell is developed. 
As we see from Figure [4~T| swimming by extending protrusions is mostly asymmetric in that they alternate 
sides, and thus the motion is not rotation-free and the trajectory of a swimming cell is snake-like rather 
than along a straight line. However we begin in Section 4.2 with a simplified symmetric amoeba, where a 


pair of side protrasions move down the cell body symmetrically so as to minimize the mechanical effects 
resulted from rotation or cell body twisting. In Section|43]we investigate how different cell and protrusion 


shapes affect the swimming ability of such a translational swimmer. Finally, in Section 4.4 we consider an 
asymmetric swimmer similar to that shown in Figure |4~T| b) and compare such a snake-like swimming style 
to the symmetric swimming style. 
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Fig. 4.1 Amoebae swim by protrusions: (a) a swimming cell with 3 protrusions jVan Haastert| 
it swims jBarry & Bretscher|2010t. 


2011 


(b) the shape of an amoeboid as 


4.1 Construction of the cell shape 


To apply the Muskhelishvili method to a swimming cell we hrst have to obtain the conformal mapping w 
corresponding to the cell shape, and then truncate its Laurent’s expansion, leaving only N negative order 
terms for some N. In general it is difficult to hnd the conformal mapping analytically that corresponds to 


a general shape, yet that of an n-polygon can be found by use of the Schwarz-Christoffel formula (Ahlfors 
[T9781 [Driscoll & Trefethen|2002| l. 

Suppose that we have an n-polygon in the z-plane with vertices zi , • • • ,Zn, and the corresponding exterior 
angles are 0i tt, • • • , 9„7r (Figure |4^. Let be the interior region bounded by the polygon. The conformal 



Fig. 4.2 The conformal mapping for the exterior region of a polygon from that of the unit disk. 


mapping from the exterior of the unit disk D in the i^-plane to is given by the following Schwarz- 
Christoffel (SC) formula ( |Ahlfors|1978||Driscoll & Trefethen|2002] l: 

z = w(C)=A + C (4.1) 

•' 9 k=i 

where Zk = w{^k) and we call l^k the prevertex to the vertex Zk under the conformal mapping z = w{Q. Since 
we are considering transformation from the exterior region of D to the exterior region of f2, the angles should 
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satisfy Y!k=\ 0k = n + 2. The prevertices can be numerically approached by using the SC toolbox given 
at http://www.math.udel.edu/~driscoll/SC/ Once we have obtained the Schwarz-Christoffel transformation 
wsc(C) for a given polygon, we truncate its Laurent expansion leaving only N negative order terms, which 
then has the form 


'^scn(C) = cti + Oo H—— I -(4.2) 

The image of the unit circle 5' under the truncated conformal mapping wscN is a contour approximating the 
original polygon. 

We approach the construction of the shape of a swimming amoeba as follows. Instead of using a polyg¬ 
onal discretization with many nodes on the boundary, we first construct an ’’inner skeleton” of the cell with 
only a few nodes (Figure |4.3| red contours), then we obtain its Schwarz-Christoffel transformation and 


Time = 0.10 Time = 0.50 Time = 0.80 



-10 0 10 -10 0 10 -10 -5 0 5 10 15 

Vertex, Time = 0.10 Vertex, Time = 0.50 Vertex, Time = 0.80 



-10 0 10 -10 0 10 -10 -5 0 5 10 15 

Prevertex, Time = 0.10 Prevertex, Time = 0.50 Prevertex, Time = 0.80 



Fig. 4.3 The conformal mapping for the exterior region of a polygon from that of the unit disk. The axes are in units of (im. 


truncate it to obtain wscn. Finally we smooth the image of the unit circle under wscn by multiplying the 
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amplitude of a\ in ( |4.2| l by a factor C S K, C > 1, to give 

+ Oo H——!-•••+ • (4-3) 

This will give us a smoothed contour enclosing the ’’inner skeleton”, and we use it as the shape of cell 
(Figure |43] black contours). Figure [43] illustrates this process, and the detailed steps of the algorithm are 
given in the appendix]^ The panels in the top row give three snapshots within one cycle. The red contours 
give the ’’inner skeletons”, each one is a prescribed polygon, with each of the four semicircle ends having 
five nodes. The black contours that are taken as the current cell shapes, are determined by a conformal 
mapping of the form of ( |4.2| i with N = 30 and C = 1.05. The panels in the middle and bottom rows show 
the distribution of vertices of the polygon (z, in Figure [4^ and the distribution of prevertices along the unit 
circle (i^, in Figure |4^ , with the correspondence relation given by the color of the dots. From Figure |43] we 
see that when the two side protrusions are close to either end of the cell body, some prevertices are crowded 
and when the side protrusions are near the middle of the cell body, the prevertices are more scattered. 

In the simulations described later we do not require strict area conservation - instead we require that the 
area changes be restricted within a small range. We define the ratio of area change within one period to be 


Ratio of area change 


Maximum of area — Minimum of area 
Average of area 


and we require that the ratio be < 0.1. 


4.2 Simulation results of swimming Dictyostelium amoebae 


We use the data for swimming amoebae from ( Van Haastert|20lT] Barry & Bretscher|2010| l. Though they 
are both Dictyostelium amoebae, they have different sizes. For simplicity we will refer to them as “van 
Haastert’s cell” (Figure [4~T](a)) and “Barry’s cell” (Figure [4~T|(b)) hereafter. Their data is shown in Table|^ 


Table 3 Experimental data for Van Haastert’s cell' Van Haastert 


2011 


and Barry’s cell Barry & Bretscher 


2010 ' 



Van Haastert 

Barry 

Maximum cell body length 

~ 25/rm 

~ 22/rm 

Average cell body width 

~ 6/rm 

~ 4/rm 

Maximum protrusion height 

~ 2/rm 

~ 4/rm 

Average protrusion width 

~ 2/rm 

~ 2/rOT 

Period of a stroke 

~ Imin 

~ 1.5min 


We use the numerical methods discussed in Section 


. 1 to generate sequences of shapes based on data 

a 


given in Table The numerical results are presented in Figure 4.4 The shapes of Barry’s cell within 
cycle for d ifferent val ues o f N are shown in Figure |4^a-c), while those of van Haastert’s cell are shown in 
Figure |4^ d-f). Figure [4^ g-i) compare the mean velocity U, mean power ^ and performance E of the two 
cells for V e [10,80]. 

A number of conclusions can be drawn from these simulations, as listed below. 
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(a) Barry, N=10 


(b) Barry, N=50 


(c) Barry, N=80 





(d) Van Haastert, NslO 


(e) Van Haastert, N=50 


(f) Van Haastert, NsSO 



Fig. 4.4 (a-c) The shapes of Barry’s cell within a cycle for different values of N. (d-f) The shapes of van Haastert’s cell within a cycle 
for different values of N. In (a-f), the top rows are the snapshots at the beginning of a cycle, the middle rows are at half of the cycle, 
and the bottom rows are at the end of the cycle. For the periods of the two cells, we have Tb = 1.5 min for BarTy’s cell and Th = I min 
for van Haastert’s cell, (g-i) A comparison of the mean velocity U, mean power 0^ and performance E of the two cells for A £ [10,80]. 


1. Protrusion shape. First we consider how the number of terms N in the conformal mappings wn affects 
the shapes of the swimmers. As we mentioned earlier, in general the term gives n angles along the 
periphery of the cell. From Figure |4.4| [a - f) we see that swimmers corresponding to wn with larger N 
have more rounded heads in the protrusions, while those corresponding to with smaller N tend to 
have sharper heads; another important difference in the shapes of the protrusions is that the connecting 
parts between the cell body and the protrusions are smoother for smaller N while more abrupt for larger 
N. Figure |43] gives an enlarged view of the protrusion regions of both cells with different A. 

2. Velocity. For either cell, the mean velocity U i ncrea ses rapidly with N when N is small, but U does not 
change much for even larger A >~ 40 (Figure 4.4 g)). For A > 40, 1/ of van Haastert’s cell fluctuates 
within a range 2.37 — 2.52fim/min while ~ 3.30/4m/min for Barry’s cell. Taking into account of our 
observation of the relation between A and protrusion shapes, our simulation results indicate that abrupt 
protrusions with rounded heads may enhance the swimming speed. 

3. Power. Unlike the mean velocity U which has a maximum as A increases, the mean power continues 
increasing as A increases (Figure |4~^h)). In particular, if we observe Figures [4~4|b,c,e,f), we see that the 
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(a) Van Haastert 



(b) Barry 



— N=10 
N=20 
N=40 
N=60 

— N=80 


pi m 


pi m 


Fig. 4.5 The areas of the protrusions in the models of van Haastert’s and Barry’s cells for different N: (a)van Haastert’s cell; (b) Barry’s 
cell. 


shapes of the same cell for N — 50 and N — SO are quite similar, yet with more terms in vv^v it requires 
much more power. 

4. Performance. Figure |4.4[ i) clearly shows that performance decreases as more terms N in W[^ are in¬ 
volved. Incorporating the observation of bump shapes, we find that smoother protrusions tend to lead to 
swimming with better performance though it might be slower. 

Comparing with experimental data, van Haastert ( |Van Haastert|201 1) reported that the swimming veloc¬ 
ity of a typical cell is ~ 3/rm/min, while our model predicts ~ 1.36 — 2.52 pm/mm. Barry at.el. (Barry 
& Bretscher||2010 l reported that the linear speed of the cells has a range of 2 — 8.4 pm/min with an 


5. 


average of about 4.2/im/min, comparing to a range of ~ 2.63 — 3.31/im/min as given by our numerical 
simulations. So far there are no experimental data regarding the power or performance of the swimming 
amoebae, and based on the data of swimming speed collected from experiments, we believe our model 
is reasonable. 

6. How cell shapes affect the swimming behavior. Both cells are the same kind of Dictyostelium amoebae, 
but from the above results and observations we clearly see that different cell shapes and sizes lead to 
the difference in their swimming behavior: as compared with van Haastert’s cell, Barry’s cell is more 
slender, with higher protrusions, and smaller in size, and the simulation results show that the average 
area for van Haastert’s cell ranges within llOprn^ — llApm^, dependin g on the value of N, while the 


average area for Barry’s cell is only within 103/rm^ — \06pm^. Figure 4.4 indicates that Barry’s cell 


exhibits faster swimming and better performance. This inspires us to study how the sizes of the cell body 
and the protrusions affect the swimming behavior of the cell, which will be discussed in detail in the 
following section. 


4.3 Effects of the protrusion height and cell body shapes on swimming 


We hrst study the effects of the protrusion shape on the swimming behavior of the Dictyostelium amoebae. 
As above we use the data for the sizes of different characteristic features of a cell in [Van Haastert] ( |201 1[ ) 
and [Barry & Bretscher (2010 1 , which are presented in Table]^ As for the conformal mappings, we truncate 
them at A = 30 so as to generate protrusions that are neither too smooth nor too abrupt. At each time step we 
adjust the cell body length to compensate for the area changes caused by the emergence and disappearance of 
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the protrusions, so as to control the area change of the whole cell within a small range. We test for different 
maximum protrusion heights, ranging from l.Sfim — Sfim. Figure 4.6 gives the relationship between the 
mean velocity, mean power, and performance and the maximum protrusion height for both Barry’s and van 
Haastert’s cells. First we see that the mean velocity increases significantly as the protrusion becomes taller 
(Figure |4~6| (a)). On the other hand the mean power for Barry’s cell increases steadily as the protrusion height 
increases, while for van Haastert’s cell the mean power first increases but for taller protrusions it approaches 
a maximum (Figure 4.6 (b)). Finally for the performance, it turns out that Barry’s cell,which is smaller, 
always performs better than van Haastert’s cell, though it swims slower than van Haastert’s (Figure |4.6| 
(c)). Moreover, while the performance for van Haastert’s cell increases as the protrusion grows higher, the 
protrusion height does not have much effect on the performance of Barry’s cell. 
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Fig. 4.6 Effects of the protrusion height on swimming for van Haastert’s and Barry’s cell: (a): mean velocity ~ maximum protrusion 
height; (b): mean power ~ maximum protrusion height; (c): performance ~ maximum protrusion height. 


Next we consider the effects of the cell body shape on the swimming behaviors of Dictyostelium amoe¬ 
bae by varying the length to width ratio of cells. We consider protrusions with the same height (^ 3/4m) 
and width (~ 2/im) and let the shape deformations have the same period (= Imin). Moreover, to reduce the 
computational effort we keep the protrusion height constant through the whole cycle, ie, we do not consider 
the emergence, growth and disappearance processes of the protrusions as in previous simulations that led 
to Figure l4^ To make the comparison fair, we control the average area of each cell within a certain range 
(~ 205 

We define the aspect ratio Rs of the cells as 

Cell body length 
* Cell body width 

so that large (small) Rs corresponds to slender (rounded) bodies. The relations of mean velocity, mean power 
and performance to the ratio of the cell body sizes Rs are given in Figure [TT] from which we see that slender 
cells swim faster than rounded ones (Figure lT^Tj ^a)), yet they require more power expenditure (Figure |4~7| b)) 
and the performance is worse (Figure |4^c)). We should note that although in general Barry’s cell is more 
slender than van Haastert’s cell, yet Barry’s cell is much more smaller, thus it swims more slowly, yet with 
better performance (Figure |43] l. 

Based on our observations on the protrusion height and the cell’s slenderness, we conclude that 

1. Within a reasonable range, protrusions with large height will result in faster swimming and better per¬ 
formance. 

2. Slender cells swim faster, but their performance is worse than those rounded ones. 
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Fig. 4.7 Effects of the cell body shape on swimming: (a): mean velocity ~ R^', (b): mean power ~ Rs ', (c): performance ~ R^. 


4.4 Do asymmetric shape deformations improve swimming? 


In Sections 4.2 - 4.3 we discussed symmetric shape deformations for amoebae swimming at LRN, but to 
date there is no clear evidence which shows that amoebae favor such symmetric modes . Rather, they seem 
to prefer asymmetric modes in which the protrusions travel one by one (Figure |4.1|l (Van Haastert||2011 


|Barry & Bretscher|2010| l. In the 2D model of swimming (Figure |44jb)) this means that the protrusions at 
the two sides alternate rather than appearing symmetrically, so that the cell swims in a snake-like trajectory. 
This raises the question as to why amoebae employ the asymmetric mode - are there advantages over the 
symmetric mode or are protrusions constrained by the internal dynamics? 

We design an asymmetric swimmer for which the protrusions alternate sides during successive cycles, 
and adapt the previous numerical scheme. Asymmetric swimmers are not rotation-free, hence we must con¬ 
sider the torque on the swimmer, and the rotational velocity is calculated via \2.21\ (Figure [47^ center). We 
compare the asymmetric swimmer with two symmetric swimmers whose cell body and protrusion shapes are 
identical, and whose conformal mapping is truncated at the same order N. One of them propagates its pair 
of protrusions at the same speed as the asymmetric swimmer (Figure |4.8| top), while the other propagates 
its pair of protrusions at half the speed of the others (Figure [4~8| bottom). The red lines in Figure |4~8]reflect 
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Fig. 4.8 A comparison of the three swimming amoebae. The top and bottom rows use the same sequence of shape changes, but 
protrusions in the bottom row travel at half the speed of those in the top row. The asymmetric swimmer in the center row alternates the 
protrusions on the two sides. The period for the symmetric swimmers is 1 min, while for the asymmetric swimmer it is 2 min. The axes 
are in units of jim. 


the trajectories of the center of mass of each swimmer, and one sees that both symmetric swimmers move 
in a straight line (i.e., rotation-free) while the asymmetric swimmer in the center images swims in a slightly 
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snake-like style. The simulations show that the symmetric swimmer propagating its pair of protrusions at 
the same speed (Figure [478] above) as the asymmetric swimmer travels a distance of X ~ 30.33 for T = 10, 
while the asymmetric swimmer travels a distance of X ^ 15.36 in the same time, which is slightly more 
than half the distance traveled by the symmetric swimmer and slightly more than that for the symmetric 
swimmer in (Figure [478] below). The average power of the asymmetric swimmer is again almost half that of 
the symmetric one with the same protrusion speed, hence their performance is the same. Thus we hnd that 
rotation that results from the asymmetric shape deformations does not lead to a reduction of performance, 
since the swimmer expends half the energy in swimming half as fast. 

Finally we show that when a protrusion always moves along one side, the cell simply rotates in the long 
run. In this case the global trajectory of the center of mass is a circle generated after a sufficient number of 
cycles (Figure [479]). 



5 Discussion 


Movement of eukaryotic cells is best-characterized for keratocytes, which use actin-driven lamellipodia 
at the leading and myosin-driven contraction at the rear to move ( Keren & Theriot|2008| l. They move with 
very little shape change, which simplifies their description, and several models that reproduce the shape have 
been proposed( |Herant & Dembo|2010|[Rubinstein et aI.|200^|Wolgemuth et a/.|201 1) |. Understanding of the 
mechanical balances that produce stable gliding motion is emerging, but the linkage between the biochemical 
state and the mechanical state is still not understood. For example, how an initially symmetric cell breaks 
the symmetry of the rest state and begins to move, and how the localization of the control molecules and the 
properties of substrate affect motion, has not been explained. 

In contrast, other cells are found to use very complicated shape changes for locomotion, and this has 
led to the overarching question posed by experimentalists, which is ‘How does deformation of the cell body 
translate into locomotion?’! Renkawitz & Sixt|2010| l. The model and results described herein on swimming 
by shape changes was motivated by recent experiments which show that both neutrophils and Dd can swim 
- in the strict sense of propelling themselves through a fluid without using any attachments - in response 
to chemotactic gradients. Our results for 2D cells show how the shape and height of a protrusion affect the 
speed and efficiency of swimming and may give insights into optimal designs of micro-robots. 

The protrusions and other shape changes used in swimming require forces that must be correctly orches¬ 
trated in space and time to produce net motion, and to understand this orchestration one must couple the 
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cellular dynamics with the dynamics of the surrounding fluid or ECM. This remains an open problem for 
future research. 
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Appendices 


A Translation of the swimmer: proof of ( |2.26| l 

Proof Uoa is defined as the translation part of the velocity field at infinity jShapere & Wilczek| 1989b) , 

f de 

Uoo = lim ® —u 

In 

where u is the fluid velocity field in the exterior domain. In 2D it can be expressed by the complex integral 


27ri f\z\^R z 


(A.l) 


Choose ^ > 0 large enough so that 0 and 'X^ are analytic on \z\ > R and continuous on |z| > R. Let ^ — z/R and then by ED and 
Table[T] 


Uoo — lim -— 




rww)-^U- 


(A.2) 


Let 


Obviously, (j >, y/ are analytic on | (^ | > 1 (including at infinity) and continuous on | (^ | > 1. Define 

= and = 

Then as functions of , 0 (1 /(^) and y/( l/(^) are analytic on | (^ | < 1 and continuous on | (^ | < 1, and on | (^ | = 1 we have 

Therefore by the Cauchy Integral Theorem 




i L. i in - «f> U = 


For the first term in |A.2| it follows by use of the Residue Theorem that 

1 r m^) 


Ini I 


=aQ. 


Moreover, since the conformal mapping z = w(^) has the form j2.15f , it is easily seen that oq = Aq and bo =Bo. Hence we have proven 
the assertion in 12.261. 
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B Rotation of the swimmer: proof of (|2.27) 


Suppose that the swimmer has the current shape and velocity field given by z = h'(^) and V((7), respectively, and denote the resulting 
torque by Tiy^w). The resulting rotational velocity 0 ) can be calculated by considering a uniform rotation of a rigid swimmer with the 
same shape z — vi^(C) resulting in the same torque. Such a uniform rotational velocity field can be expressed as 

y((7) = icow{a) 


Let 


y™^(cT;r) ^ /w(fT;r) 

i.e., V — coV^°^. To match the torque that results from the two velocity fields, we have T(V;w) — Thus co can be expressed 

as 

T{V;w) 




which gives the first relation in \2.21\. 

For the other two relations in \2.21\, first we show that b-\ —B-ya-i as follows. From (|2.13t it is easily seen that 


b-\^ 


2ni Jd £2 

Since dQ — {w((7;r); <7 G 5^}, the above equation can be transformed into 


^ / Viz)‘iz. 

J dQ 

iformed into 


(B.l) 


\dG — 


Finally from equations i |2.15|2.18|.B.l) we have 

, 1 f , B-i B-2 n/ Ot-i 2a-2 

2m 7si ^ a (j2 ^ ^ <7^ (7^ 

Finally we prove the relation T{V;w) — —47ifi3b_i as follows. 

Proof The torque associated to the boundary condition V{<7) is given by 

T{V\w) — lim 3 ^r X fds — lim 3 J 

From Table^we see that / is a sum of two parts: along the n direction and —along the « direction, where 

n — —idz/ds is the exterior normal on dQ. When taking the cross product with r, the first part necessarily vanishes since it is parallel 
to r. Hence 


and 


On the other hand we have 


^rx fds ^ -2i^ z(f:(l>" + dz 

r(y;w) ^- 2/1 lim 9^ / z(z(!)" + X^^^dz 

R^<=o J\z\=R \ / 


|“ = -(zr+x') 


r(y;w) = 2^ lim / z^dz = 2)1 lim — f udzl 
J\z\=R dz L J\z\=R J 

= lim 3 [ T —: [ u{z,z)d^ . 

R — \-27tl J 


SO 









27 


Now we only need to calculate the complex integral in the above equation. We proceed as in the proof for For R large enough, 
we introduce the substitution ^ =z/R and then we have 


T{V-w) 

4kii 

On 1I = 1, we have so 




T(V-,w) 

47111 


= fc3[- WlL, (^- M)- pVM))^4 


By the Residue Theorem 




27liJ\^\=l ^ 


2 ^^=0 


and by the Cauchy Integral Theorem 


i L. f = i in, =J'f'U. ■ =“ 

where 0^(oo) — 0 because of the form of (j) given in Thus the first two terms vanish, and by the Cauchy Integral Theorem, 

i .in, i i m^] Uo= i k(|)] Uo 


By j2.13k, on 1^ I < 1 we have 


thus 


Finally we have that 


Z'( = )=5o+'^^ + ^^^+ - + '^«" + - 






T{V-,w) 

471)1 


= 3 R^ =-3h- 


C Algorithm for the shape changes 

The inner skeleton, i.e. the red contour in Fig 4.3, is a polygon of cross-like shape, with a semi-circle at each end of an arm. The polygon 
in each step consists of 28 vertices, where each semi-circle end has 6 vertices. The initial shape pO (i.e, shape at t=0) is generated as 
follows. 

p0(l)=2-l-2i 

p0(7)=17.5-2i 

p0(8)=17.5-2.2i 

p0(14)=19.5-2i 

p0(15)=20-2i 

pO(21)=I9.5-i-2i 

p0(22)=19.5-i-2.2i 

pO(28)=I7.5-i-2i 


for k=l:l:5 
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pO(l+k) = 2+ 2i*exp(k*li*pi/5) 
pO(8+k)=18.5-2.2i - l*exp(k*li*pi/5) 
pO( 15+k)=20-2i*exp(k* 1 i*pi/5) 
p0(22+k) = 18.5 + 2.2i + l*exp(k*li*pi/5) 


end 

pO(29)=pO(l) 

For each time step dt, move the two arms (i.e., the blebs) of the cross forward simultaneously. We change the length of the blebs 
and the body so as to simulate the grow and decay of the blebs, and compensate for the resulting area change. Below is the pseudo-code 
that generates the inner skeleton p at the k'* time step. 

for j=l:l:6 

P(j) = pO(j) + 0.5*H*sin(dt*k*pi) 

end 

p(7)= p(7) - 0.5 
for j=8:1:13 

P(j) = p0(j)- 0.5*k - H*sin(dt*k*pi)*li 

end 

p(14) = p(14) - 0.5 
forj=15:l:20 

P(j) = p0(j) - 0.5*H*sin(dt*k*pi) 

end 

p(21) = p(21)-0.5 
for j=22:1:27 

P(j) = p0(j) - 0.5*k + H*sin(dt*k*pi)*li 

end 

p(28) = p(28) - 0.5 
P(29)=p(l) 
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